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Abstract - 

A model of rods with heads of variable size, which are confined to a planar 
surface, is used to study the influence of the head group size on tilted phases 
in Langmuir monolayers. Simple free energy considerations as well as exact 
zero temperature calculations indicate that molecules with small head groups 
tilt towards next nearest neighbors, and molecules with larger head groups 
towards nearest neighbors. This provides a possible explanation for recent 
experimental results, and for details of the generic phase diagram for fatty 
acid monolayers. 



PACS numbers: 68.18, 68.35 



1 Introduction. 



Monolayers of simple amphiphiles at the air-water interface form a complex 
variety of condensed phases |], ||. A generic phase diagram for fatty acids is 
shown in Figure 1. The phase behavior appears to be driven mainly by the 
conformations of the aliphatic tails of the molecules. However, head group in- 
teractions also play a role, especially in tilted phases. Shih et al @ and Fischer 
et al M have recently performed systematic studies of mixtures of heneicosanic 
acid and heneicosanol. They find that the L 2 (h,d) phases with tilt towards near- 
est neighbors are gradually suppressed, as the alcohols are added to the fatty 
acid monolayer, until they finally disappear at alcohol concentrations above 
35 %. In pure heneicosanol monolayers, only phases with tilt towards next 
nearest neighbors exist. 

According to a common picture, tilt order in monolayers is often induced 
by a mismatch between the head group spacing and the chain diameter. The 
size of the head groups plays a crucial role for this mechanism. From the 
observation that the zero pressure lattice constant is smaller in heneicosanol 
monolayers than in heneicosanic acid monolayers |§, one can deduce that the 
effective size of alcohol head groups in water is smaller than that of fatty acid 
head groups. The results of Shih et al and Fischer et al therefore suggest that 
the head group size also determines the direction of tilt. Further investigations 
of surface pressure - temperature phase diagrams for monolayers of fatty acids 
with esters by Teer et al || are consistent with this picture. 

We investigate this possibility in the framework of the simplest possible 
theoretical model: Amphiphilic molecules are represented by rigid rods with 
heads larger than the rod diameter. We discuss the model first using a simple, 
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general free energy ansatz. Our results are confirmed by exact zero temperature 
calculations for a system with Lennard- Jones interactions. 

2 Free Energy Considerations. 

We consider rigid rods of length L with steric repulsive interactions and longer 
range attractive interactions; the favorite distance between two rods at given 
temperature is r t . The rods have heads of radius r h , which are confined to 
a planar surface and arranged on a distorted hexagonal lattice, with variable 
lattice constants a and b (Figure 2). Two cases are discussed, tilt in the 
direction of a (nearest neighbors, NN) and tilt in the direction of b (next 
nearest neighbors, NNN). In this geometry, rods have two different types of 
neighbors, labelled (i) and (ii) in the following. 

The head size gives rise to the geometrical constraints 

(i) a>r h (ii) (a 2 + b 2 )/A>r 2 h . (1) 

If the heads are small, rh < r t , the optimal configuration is one of untilted 
rods arranged on an undistorted hexagonal lattice, with the lattice constant 
a = r t . The volume per rod of the rod layer is then given by V = y/3/2 r\ L. 
If the heads are larger than the tails, the attractive interactions of the tails 
cause them to tilt. This may involve a distortion of the lattice; the surface per 
rod A occupied by the monolayer increases, which is in turn penalized by the 
surface tension E and the surface pressure II. 

Such considerations lead to the free energy ansatz 

F = vj-(^-ir + (j: + u)A (2) 

with the volume compressibility of the rods k and the volume per rod V. The 
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repulsive interactions between the rods cause additional constraints 

a^r t (a 2 + & 2 cos 2 #)/4^r 2 (NNN) 

W acos#~r t W (a 2 cos 2 fl + 6 2 )/4^r 2 (NN) U 

where # is the tilt angle. It is convenient to define the mismatch parameter 
5 = 1 — Th/rt, and to rewrite the free energy as 

F oc F = A(M- X [Z) 2 + — with A = — ^ r. (4) 

Our task is to minimize F with respect to the surface A = ab/2 and volume 
V = AL cos 8 > Vb, under the constraints (0) and (^|). The remaining indepen- 
dent model parameters are A and 5. Assuming that k = A5 is of order unity 
or less, the minimum value of F can be expanded in powers of 5. Details of 
the calculation are presented in the appendix. 

In layers with rods tilted towards next nearest neighbors, the volume per 
rod cannot be optimized, because the distance to the neighbors (i) is too large, 
a > Th > r t . The free energy is minimized at A = \pir 2 h and takes the value 

F NNN = 1 + (2 + h)6 + (1 - ^k)6 2 + ■■■ (5) 

The head lattice is not distorted. This is an effect of the "hard" constraints 
of eqn (|I]), i.e. the hard core interactions between heads. If the heads interact 
via a softer potential, the lattice gets slightly distorted in the direction of the 
tilt. 

In layers with rods tilted towards nearest neighbors, the volume per rod 
can be optimized, but at the cost of an increased surface energy. Surface and 
volume contributions have to be balanced. For large k, k ^> 1, the volume 
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term determines most of the behavior, and one gets 

p fc»i / 1 + 4:5 - 65 2 + • • • : 8 > 5 (a) 

™ ~ \ 1 + (4 - 1/Jfe)<5 - (6 - 16/Jfe + 10/fc 2 )5 2 + • • • : 5 < 5 (b) 

(6) 

with 5q = \JS/2 — 1 = 0.225. In this case, the head lattice is distorted, and the 
volume per rod is minimal (V = Vq) in the limit k — > oo or 5 > Sq. For small 
on the other hand, fc< 1, the surface term dominates. The heads form an 
undistorted lattice, and the free energy is given by 

F NN k ^ l + (2 + k)6+(l- k)5 2 + ■■■ (7) 

In the intermediate regime, k ~ 1, the smaller of the two solutions (||) or (0) 
applies. Note that the free energy (^) is always larger than ([5]). Hence rods 
tilt towards next nearest neighbors at k <C 1, or generally if distortions of the 
head lattice are not allowed. 

Suppose now that A is kept fixed, and the head size is gradually increased. 
Rods with small heads, 5 < 8 c i = 0, are untilted. If the heads are just slightly 
larger than the tails, < 5 < 1/A, the case k < 1 applies, and the rods tilt 
towards their next nearest neighbors. On increasing the head size (and if A is 
sufficiently large), a swiveling transition to a phase with tilt towards nearest 
neighbors takes place. The transition point <5 c2 is located by equating eqn (^) 
with eqn @. One thus gets a sequence of phase transitions at 

5 cl = and 5 c2 = k c /A with k c = 3.93 (8) 

to lowest order of 5 (A > k c /8o). Note that the transition at 5 C \ from the 
untilted to the tilted state, does not depend on A. This is again a consequence 
of the hard core interactions between the heads. If the heads are soft, they are 
squeezed together at high pressures {i.e., small A), and the rods stand up. Note 

5 



also that the free energy (0) is temperature dependent via the temperature 



dependence of the compressibility k and the surface tension E. However, the 
model cannot be expected to reproduce thermal properties of a monolayer, 
since the hard core interactions between the heads are basically athermal. 

3 Exact Ground State Calculations. 

In order to test these predictions, we have studied numerically the ground 
state of a system of rods with Lennard- Jones interactions between tails and 
repulsive soft core interactions between heads. The interaction energy of two 
rods at grafting distance r, with tilt direction e L , is given by 



where d is the distance between elements dl and dl' on the rods, d = \ f+ (I — l')eL 
and the potentials are taken to be 



The potentials Vlj and Vh are cut at d = 2 and r = a, respectively, in order 
to make the problem more tractable. The total energy E t is the sum over all 
nonzero pair interactions, i.e., the interactions of a rod with up to 24 pairs 
||. It is minimized with fixed molecular area A by numerical methods. The 
surface pressure is obtained using II = —dE t , m inl dA. 

The lattice constant and the compressibility of a system of infinitely long 
rods can be calculated analytically. We identify these parameters with the rod 
diameter r t and the volume compressibility k (r t = 0.9333 and k = 0,0184), 




(9) 




(10) 



(11) 
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and a with the head size rv This allows a quantitative comparison of the 
numerical results with the predictions of the previous section. 

The phase diagram for chain length L = 5 is shown in Figure 3. With 
increasing head size, one finds the sequence of phases predicted in section 1: 
Untilted (U), Tilt towards next nearest neighbors (NNN) and towards near- 
est neighbors (NN). Due to the soft core interactions between the heads, the 
first transition U-NNN becomes pressure dependent. Apart from this detail, 
the results are in qualitative and reasonable quantitative agreement with the 
predictions of eqn flj|). 

4 Discussion. 

In sum, we find that tilt order can be induced in systems of rods with grafted 
heads by a mismatch of head size and rod diameter. Rods with large heads tilt 
towards nearest neighbors, and rods with heads just slightly larger than the rod 
diameter tilt towards next nearest neighbors. The critical head size, at which 
the swiveling transition occurs, increases with pressure and decreases with the 
rod length. If distortions of the hexagonal head lattice are not allowed, rods 
always tilt towards their next nearest neighbors. 

This is in agreement with previous theoretical predictions. Scheringer et al 
report that rods grafted on an undistorted hexagonal lattice tilt towards 
next nearest neighbors. Kaganer et al || find that this remains true, even if 
one allows for lattice distortions, in systems of pure rods (no heads) which are 
forced to tilt by the constraint of a fixed homogeneous surface density (see also 



ID |). Some authors study the effect that tilt towards nearest neighbors can be 



induced by additional, attractive interactions between the rods and the surface 
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H [|. The present work indicates that it might simply be a consequence of 
the larger head group size. 

Our effect offers a simple explanation for the observations of Shih et al || 
and Fischer et al discussed in the introduction. It also provides an inter- 
pretation for earlier experimental findings. In particular, the fact that phases 
with tilt towards next nearest neighbors appear at high surface pressures (see 
Figure 1), comes out quite naturally. 

Shih et al [IT] have studied monolayers of heneicosanic acid in the presence 
of calcium ions in the subphase, and systematically varied the pH of the sub- 
phase. At low pH, they find a sequence of phases with tilt towards nearest 
neighbors, tilt towards next nearest neighbors, and no tilt, as a function of 
pressure. On increasing the pH, the phase tilted towards next nearest neigh- 
bors moves down to lower pressures, and finally replaces the phase with tilt 
towards nearest neighbors. At even higher pH, the monolayer is untilted at 
all pressures. Assuming that the effective size of the COOH - head groups is 
reduced at high pH - e.g., due to more efficient screening of the electrostatic 
interactions at higher concentration of positive ions in the subphase - this 
result corresponds exactly to the predictions of our model. 

To conclude, our simple considerations provide an explanation for a remark- 
able variety of experimental observations. Obviously, the real interactions be- 
tween head groups in Langmuir monolayers are much more complicated than 
those assumed here. Nevertheless, the success of our model suggests that the 
influence of head groups on the microscopic structure of Langmuir monolayers 
is to a large extent simply determined by their size. 
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Appendix 

We wish to minimize the free energy @ with the constraints (|l|) and 
For tilt towards next nearest neighbors, the solution can be written down 
immediately: a = r h , b = V3r h , cosd = yj (A — (1 + <5)) 2 )/3(l + 5) 2 , and 

^iVAW = (1 + o J + A( — , = - 2 ) i 12 ) 



V3(1 + 5)V4- (1 + 5) 2 
The expansion of F NNN in powers of 5 yields (|5|). 

For tilt towards nearest neighbors, different solutions are possible. If the 
volume term dominates, one achieves V — Vo with cos# = r t /a, b = V3r t and 



minimizes the surface A under this constraint with a/r t = y 4(1 + 5) 2 — 3. 

F£2 = <M + ^) 2 " 3. (13) 



In this case, the head lattice is distorted. Expansion of F$2 gives (f|a). 

If the surface term dominates, the head lattice is not distorted, (6 = Vor^, 
a = and with this constraint the volume is minimal for cos 8 = r t /a. 

4 A i=(l + ^ + A^y (14) 

Expansion yields (|7|). 

If none of the two terms dominates, the solution lies between (|13|) and (|T~ 



b = y/3r t y and a/r t = dA(l + 5) 2 — 3y 2 with 1 < y < 1 + 5 (cos 8 = r t /a). 



F NN = min F(y) with F(y) = yd Ail + 5) 2 - 3y 2 + '- (15) 

Ky<l+8 v y 



(y-i) 
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Note that dJ- / dy\ y =i^S) for 5^5$ = \Js/2 — 1. If 5 > 5q, 5q, the function 



J-{y) has no minimum in the interval ]1, 1 + 8[, and Fnn is either given by 
= Pjfy, or JF(1 + 5) = i^/v- If 5 < 8 , J-(y) has a minimum. The surface 
term still dominates at small A (A ~ 1/5), F/viV = Fjfy. At larger A, no closed 
expression for F^n can be given, but an expansion in powers of 5 is possible: 
The minimum of T{y) is found at y = 1 + a8 with 
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21 - 16Jfe 
2A; 2 



5 + 



4(40 - 42A; + 9k 2 




(16) 



(fc = A5), which leads to eqn (|.b). 
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Figure Captions 

Figure 1: Part of the generic phase diagram for fatty acid monolayers. Solid 
lines are first order transitions, broken lines are second order. The am- 
phiphile molecules are tilted towards next nearest neighbors in the phases 
Ov and L' 2 , towards nearest neighbors in L 2 d and L 2 h, and untilted in S 
and LS. In the S, L 2 and L^h phases, the backbones of the chains are 
ordered, the other phases are rotator phases. Increasing the length of the 
aliphatic chains shifts the phase diagram to higher temperatures (after 

0)- 

Figure 2: Side view of the rigid rod model (top) and top view of the head 
lattice (bottom). We consider tilt directions in the direction of a (nearest 
neighbors), or b (next nearest neighbors, (i) and (ii) mark the different 
types of direct neighbors in this geometry. 

Figure 3: Ground state phase diagram of the model of section 3 (rods with 
Lennard- Jones interactions) in the plane of head size a vs pressure II, at 
rod length L = 5. Phases are: tilt towards next nearest neighbors (NNN), 
towards nearest neighbors (NN), untilted (U). Dashed lines indicate the 
prediction of eqn (^) with E = 6.7. 
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